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BERNOULLI NUMBERS AND SOLITONS 


M-P. GROSSET AND A.P. VESELOV 


Abstract. We present a new formula for the Bernoulli numbers as the fol¬ 
lowing integral 


B2m — 


(- 1 )" 


2‘2m+l 


/ -l-oo —1 

(- -sech^x)^dx. 

'da;—1 


This formula is motivated by the results of Fairlie and Veselov, who discovered 
the relation of Bernoulli polynomials with soliton theory. 


Dedicated to Hermann Flaschka on his 60*^ birthday 


1. Introduction 

In the paper ^ D. Fairlie and one of the authors discovered an interesting 
relation of the Bernoulli polynomials with the theory of the Korteweg-de Vries 
(KdV) equation 

Ut - 6uUa: + Uxxx = 0. 

It is known since 1967 due to Gardner, Green, Kruskal and Miura |2] that this 
equation has infinitely many conservation laws of the form 


I [tt] Dm (t 


T)dx, 


where Pm are some polynomials of the function u and its ^-derivatives up to order 
m. They are uniquely defined by some homogeneity property modulo adding a 
total derivative and multiplication by a constant. This constant can be fixed by 
demanding that Pm{u, Ux, Uxx, Um) = Um a function of derivatives of order 
less than m. 

The KdV equation is famous for its remarkable family of solutions, known as 
solitons, the simplest of which is a one-soliton solution 

u = —2sech^(a; — 4t), 


corresponding to the initial profile u(x, 0) = —2sech^a;. 

The main result of PP is the following formula relating the Faulhaber polynomials 
Fm with the integrals of the KdV equation: 


( 1 ) 


/m-i [—2Asech^a:] = (—1) 


Fm{\). 


2m -I- 1 


Recall that the Faulhaber polynomials are directly related to the Bernoulli polyno¬ 
mials through the formula 


( 2 ) 


B 2 m+ 2 {x -|- 1) — {2m 2)Fm 
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where Bk{x) and Bk are Bernoulli polynomials and Bernoulli numbers respectively 
(see 012 ) The Bernoulli numbers have the following generating function: 


z 


- 1 


E 




k\ 


All odd Bernoulli numbers except Bi = —^ are zero and the first even Bernoulli 
numbers are 


Bq = 1, B 2 


1 

6 ’ 


B 4 



42’ 


B. 


_1 B =1 B 691 

30’ 66’ 2730 


They play an important role in analysis, number theory, algebraic topology and 
many other areas of mathematics. 

In this note we show that the relationship with soliton theory brings the following 
formula for the Bernoulli numbers: 

/_i'\m —1 /•+00 

(3) B 2 m = ^ J ((sech^a;)(™“^))Va;, m > 1, 

where (sech^a;)^’”"^^ denotes the (m — l)th derivative of sech^x. We present also a 
direct proof of this formula based on the elegant arguments due to Logan, which 
we found in the book by Graham, Knuth and Patashnik. 


2. Proofs of the main formula. 


Let us first present the proof based on the results of Q] . Recall that the Faulhaber 
polynomial Am (A) with m > 1 has the form 

IP ^ \ 2 I ^m\3 I I \m+l 

Fm(A) — Q!2 A + a3 A + ... + 

with some rational coefficients a™, a™,..., (see 13). Since the only quadratic 
term in the density of the KdV integral Im-i is / u^_idx, we have from the relation 
that 


( 4 ) 


/ + CX3 

-00 


2m + 1 


Differentiating the formula 0 twice with respect to x gives 


^ 2x+1 ^2t^„ + ( X{x + l) ^ _ _1 

\ o ^m\ t-> )'^m\ o o I 


yB 2 m +2 (l+x) — (2m+l)i32m(l+2^) 


2 ' 2 ' 2 ^ 2m + 2 

because B 2 m+ 2 {x)" = (2m + 2)(2m + \)B 2 m{x). Since F^iO) = 0 this reduces to 

iA"(0) = ( 2 m+l)R 2 m(l). 

Now using the well-known symmetry Bk{l — x) = (—l)^i?fc(x), we have i? 2 m(l) = 
A 2 m( 0 ) = B 2 m and thus 


a™ = ^^’m(0) = 2(2m+l)R2m. 

Substituting this into 0 we come to the formula ©■ 

We are now going to prove the formula 0 directly without reference to soliton 
theory. We borrow the main idea from the book jS], where it is attributed to Logan. 
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Consider the integral Jm = ^ ^^)‘^dx. Integrating by 

parts m — 1 times gives 

1 ^ + oo 1 ^ + oo 

Jm = 22 m+i J (sech^a:)^^'"~^^secUa:(ja: = 22771+1 J tanhcc^^^^^Hanhx^^Ux. 

Let y = tanh x then 


J m. 


1 


r +1 


22m+l 


T2m-i{y)dy, 


'-1 


where the polynomial Tk{y) is the fc-th derivative oiy = tanh* rewritten in terms 
of y: 

Ti = = 1 - 

T 2 = = - 2 yy(i) = - 2 y(l - y^) = - 2 y + 2y^ 

Ts = y^^^ = -2y^^'> + 6y^y(^) = (1 - y^)(-2 + 6y^) = -2 + 8y^ - 6y^ 

T 4 = y^^^ = 16yy^^^ — 24y^y^^^ = (1 ~ 2/^)(16y — 24y^) = 16y — 40y^ + 24y® 
and so on. 

These polynomials can be determined by the recurrence formula 
(5) T„(*) = (1-*2 )t„_i(*)' 

with To(*) = X. ^ They have integer coefficients with the highest one being equal to 
(—l)"n! and have the symmetry T„(—*) = (—1)"“^T„(*). Note that the non-zero 
coefficients have the alternating signs and their total sum is zero. The integers 
T' 2 m-i( 0 )| are called tangent numbers (see e.g. 0 ). 

Let us consider now the generating function 

Tix,z) = '^Tn{x) — . 

n\ 

n>0 


Lemma 1. The generating function for the polynomials Tn{x) is 

, sinh z-I-* cosh z 

( 0 ) T{x,z) = -^-——. 

cosh z + X smh z 

Indeed it is easy to see that when x = tanhw the function T{x,z) becomes 
tanh(z -|- w). Now the claim follows from the Taylor formula and the definition of 
the polynomials T„(*). 

From this one can derive an interesting relation between Bernoulli and tangent 
numbers, see [H], formula (6.93) and discussion after that. We will need however a 
slightly different result. 

Lemma 2. Bernoulli numbers can be written as 


(7) 


Bm = 


I 


2 m+l 


Tm-i{x)dx 


for all m > I. 

Indeed, this formula is obvious for odd m bigger than 1 since Tm-i{x) is an odd 
function and both sides of o are then equal to 0. In order to prove the formula 


^These polynomials should not be confused with the classical Chebyshev (Tchebycheff) poly¬ 
nomials, which are also denoted as Tn{x). We follow the notations of the book |2|, where trigono¬ 
metric version of these polynomials was considered. 
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for even m, let us consider first the left hand side of The function T{x,z) can 
be rewritten as 


sinh z 


, , coshzlcc-i-r- 

T{x, z) = - ^ 


cosh 2 > 


sinhzlcc ■ i. , 

^ sinh z > 

Integrating with respect to x gives 


= coth z — 


1 


sinh z(x + cothz) 


which leads to 

( 8 ) 


J T{x, z)dx = X coth z— 
j T{x,z)dx = 2cothz— 


sinh^ 


■ In I 


coth; 


sinh^ z 


= 2 coth z + 2z coth^ z. 


The expansion of the function coth z can be written in terms of the Bernoulli 
numbers: 

o2n TD 


(9) 


1 


1 ^ \^ 

cothz =- h 2_^ 0,2n-lZ 


2n-l 


) a2n-l — 


n—1 


(2n)! 


This easily follows from the identity 


z z z z 

— coth — = — I-. 

2 2 2 - 1 


From ® and m it follows that 


( 10 ) 


J T{x, z)dx = 4 na 2 n-iz^ 


On the other hand from we have 

(11) f T{x,z)dx = '^{ f T 2 n-i{x)dx) 


n>l 


^2n-l 


(2n- 1)!’ 


since T 2 n(x) are odd polynomials. Comparing m and CD and replacing a 2 n-i 
by we obtain 


Bon — 


1 

2^ y_^^2n-i(a^)dx. 


This proves Lemma 2. 

Remark. As a corollary we have the following representation of the Bernoulli 
number Ran as a sum of the fractions of the type 


Ran — 


22ri 


h + h + h 

1 3 5 


krt 


2tl 1 


where the denominators are consecutive odd numbers and the numerators are the 
coefficients of the polynomials T 2 „_i(x). Comparison with the formula (6.93) from 
the book [5] leads to the following relation with the tangent numbers: 


( 12 ) 


j ^ Tn{x)dx = Tn(0). 


2 "+! _ 

Now combining Lemma 2 with the definition of the polynomials we have our 


mam 
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Theorem. 

representation 


Bernoulli numbers B 2 m with m > 1 have the following integral 

^—1 /•+00 ^] m—l 


^2m — 


22 m+l 


/ +CX) ^m—L 

(- - Tsech^xYdx. 

-oo rfai—1 ’ 


3. Concluding remark. 


There exists a generalisation of the Bernoulli numbers, known as Bernoulli- 
Hurwitz numbers BHk and related to the coefficients of the Laurent series of the 
Weierstrass elliptic function p{z) at zero: 


(13) 


piz) = z ^ + 


BH2k-\-2 2k 

{2k + 2){2k)\ 


(see inilZl)- 

Our formula suggests a different generalisation of Bernoulli numbers related to 
an elliptic curve, namely 

(14) Bfl = I (p(z)(™-i))2dz, m > 1, 


where the integral is taking over one of the cycles on the corresponding elliptic 
curve. Obviously these numbers depend also on the choice of such a cycle, which 
is an element of the first homology group of the curve. They naturally appear in 
relation with the spectral density for the classical Lame operator [^. 
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